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The cosmological power spectrum of the coherent matter flow is measured exploiting an improved 
prescription for the apparent anisotropic clustering pattern in redshift space. New statistical analysis 
is presented to provide an optimal observational platform to link the improved redshift distortion 
theoretical model to future real datasets. The statistical power as well as robustness of our method 
are tested against 60 realizations of 8 Gpc^ dark matter simulation maps mocking the precision 
level of upcoming wide-deep surveys. We showed that we can accurately extract the velocity power 
spectrum up to quasi linear scales of fc ~ 0.1 /iMpc~^ at z — 0.35 and up to ~ 0.15 /iMpc~^ at 
higher redshifts within a couple of percentage precision level. Our understanding of redshift space 
distortion is proved to be appropriate for precision cosmology, and our statistical method will guide 
us to righteous path to meet the real world. 
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I. INTRODUCTION 



The emergence of a standard model for the Universe 
dominated by an unknown substance of dark materials 
has revolutionized our understanding of the Universe. 
Since the first firm evidence of dark energy in 1998 [H, [3, 
there has been substantial observational and theoreti- 
cal research aiming at understanding the true nature of 
this phenomenon. In recent years, many authors have 
started exploring the possibility that dark energy, and 
the observed acceleration of the expansion of the Uni- 
verse, could be the consequence of an incomplete theory 
of gravity on cosmological scales and may require modi- 
fications to Einstein's theory of General Relativity. The 
information of underlying science about the Universe is 
given by looking at structure formation on large scales. 

Several authors have shown that by combining various 
probes of the large-scale structure in the Universe, it is 
possible to test the relationship between these quantities 
which, in the linear regime, can generally be described by 
two functions of time and scale . Those can be con- 
strained through a combination of weak Icnsing and an 
independent probe of matter-energy fluctuations. This is 
motivated by the fact that the weak lensing experiments 
probe the geometrical potential of combination between 
curvature perturbation and Newtonian force, which de- 
termines the trajectories of photons through the Uni- 
verse, while matter fluctuation measurements probe the 
Newtonian force alone determining local inhomogeneitics 
of matter-energy. 

The coherent motion of the galaxies opens a unique 
opportunity to access the fluctuations of the underlying 
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matter density field [9l-l24l|. This technique relies on the 
redshift space distortions seen in galaxy surveys. Even 
though we expect the clustering of galaxies in real space 
to have no preferred direction, galaxy maps produced by 
estimating distances from redshifts obtained in spectro- 
scopic surveys reveal an anisotropic galaxy distribution. 
The anisotropies arise because galaxy recession velocities, 
from which distances are inferred, include components 
from both the Hubble flow and peculiar velocities driven 
by the clustering of matter |25l - [27| . Measurements of the 
anisotropies allow constraints to be placed on the rate of 
growth of clustering. 

Measurements of coherent motion field from redshift 
distortion maps have been plagued by systematic uncer- 
tainties which have made their cosmological constraints 
uncompetitive compared to other probes of the Universe. 
The cosmological density and velocity field couple to- 
gether and evolve nonlinearly. In addition, the mapping 
formula between the real and redshift space is intrinsi- 
cally nonlinear. These nonlinearities prevent us from 
inferring the linear coherent motion from the redshift 
space clustering straightforwardly. Recently, an accu- 
rate theoretical model for the redshift distortion was pro- 
posed in [2^. They take into account the fact that the 
linear squeezing and non-linear smearing effects on dis- 
torted maps are not separable to each other and develop 
a more elaborate description than simple factorized for- 
mulation. The derived correction terms at leading higher 
orders assist us to achieve better fit to simulated data. 
They also include non-linear corrections formulated us- 
ing the closure approximation to predict the non-linear 
growth in density-density, density-velocity and velocity- 
velocity spectra. In this work, we assume a perfect cross- 
correlation between density and velocity fields at linear 
level to decompose the coherent motion spectra. Thus 
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sity power spectrum and its velocity counterpart. 

Despite all theoretical efforts, the exact form of FoG ef- 
fect is unknown. We adopt most common functional form 
of FoG effect, such as Gaussian or Lorentzian, and pa- 
rameterize FoG effect using one-dimensional velocity dis- 
persion representing the randomness of the motion which 
erases the correlation structure on small scales. The pa- 
rameter space is extended to include this uncertainty of 
FoG effect in addition to the linear spectra of density and 
velocity fields, which were originally proposed by [29l.[30j. 
We run MCMC routine to find best set of spectra and 
FoG parameter, and we find that coherent motion spec- 
tra are measurable at linear regime with good precision. 
Coherent motion spectra are measured at a couple of per- 
centage accuracy with appropriate k cut-off. 

The paper is organized as follows. In Sec. |TT1 we begin 
by introducing suggested theoretical models of redshift 
distortions. We then describe the statistical method to 
extract the coherent motion spectra. Using this method. 
Sec. mil presents our main results on the measurements 
of linear density and velocity spectra. In Sec. IIVI the 
impact of wrong prior assumption on the decomposition 
of power spectra is discussed. Finally, we conclude in 
Sec. El 



II. METHODOLOGY 

We first highlight analytical models for the redshift- 
space power spectrum. We then show our methodology 
of reconstructing the linear density and velocity power 
spectra from the two-dimensional power spectrum ob- 
served in redshift space. 



Analytical models for the power spectrum in 
redshift space 



At large scale, we expect that the density field as well 
as the velocity field are small perturbations to the homo- 
geneous universe. When the higher-order contributions 
are negligibly small, the two-dimensional power spectrum 
in redshift space, P{k,fi), can be expressed as 

P{k,f,) = + 2A.2p]g(fc) + ^^^Pteik), (1) 

where we denote by 6 and O the density contrast and 
the velocity divergence, with the latter defined by 6 = 
— (1-1- z)\/v/H. The auto- and cross-power spectra of 
the two fields in linear theory are expressed as P}j^{k), 
with i and j being either S or Q. This formula describes 
the effect of coherent velocity fiow towards overdensity 
at large scale (Kaiser effect). Because of this effect, the 
clustering pattern in redshift space is enhanced along the 
line-of-sight direction. If the gravitational law follows the 
general relativity, the three spectra are not independent. 



but related as 
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where f{z) = dlnZ?+(z)/dlna is the growth rate param- 
eter with being the linear growth rate. In this case, 
all the velocity information is contained in the parameter, 
/(z), which is constant over wavenumber k. 

However, the planned/ongoing galaxy redshift surveys 
as well as the existing large surveys mainly target weakly 
nonlinear scale, where the feature of BAOs is prominent. 
Moreover, by appropriately modeling this regime, we can 
in principle enlarge the range of wavenumber to be taken 
into account in the analysis, and improve the constraints 
on the gravitational law. Thus, we have to somehow 
incorporate nonlinearity to make maximum use of these 
surveys. 

First of all, the cross- and auto-power spectra of the 
density and the velocity fields are naturally expected to 
receive nonlinear corrections. Another important effect 
arises from the random motion of galaxies. By combining 
these two effects, Scoccimarro [l^] proposes the following 
formula 

P{k,^l) = {P55(fc) + 2/i2p,-e(fc) + /Pee(fc)} G^(M,(4) 

where the effect of random motion is captured by the fac- 
tor G^{kn). Note that we have replaced the three linear 
spectra in Eq. ([1]), Pjj^^k), with their nonlinear counter- 
parts, Pij{k). In the original paper by Scoccimarro [27| . 
the factor G'^(fc) is designed so that it accounts for the 
random motion of the galaxies at large scale, and he con- 
siders the Gaussian shape for this factor: 



G^ikfi) = G^^^ikfi) = cxp {-(fc^cr^)2} , 



(5) 



where (Tv denotes the dispersion of the one-point PDF 
of the velocity in one-dimension. At smaller scale, in- 
side the cluster of galaxies, the virial motion of galaxies 
also give a suppression of the power spectrum in redshift 
space. This effect is called Finger-of-God (FoG), and can 
also be approximately described by multiplication of the 
factor G^{k(x). A Lorentzian form of this factor has been 
frequently adopted based on the results of A^-body sim- 
ulations: 



G^{k^l) = G^"«(fc/.t) 



1 + {k^ia^f 



(6) 



For convenience, we do not distinguish between these 
damping effects, and simply call them as FoG in this pa- 
per, although, strictly speaking, the former one at large 
scale has a different origin. 

Recently, Taruya, Nishimichi & Saito [2^ proposed a 
more accurate model for the redshift-space distortion. 
Motivated by the fact that the Kaiser and the FoG effects 
can not be separated, and should not be described as a 
factorisablc form as in Eq. (U), they derived correction 
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terms which have long been overlooked. Their formula 
reads 

+A(fc,^)+i3(fc,Ai)}GX(M, (7) 

where the full expressions for the terms A{k,fi) and 
B(k,fi) based on perturbation theory can be found in 
the Appendix A of that paper [2^. Let us note some 
important features in the new terms. First, they include 
higher-order polynomials in fi and /: the A-term has a 
term which scales as f^^J,^, while we have a /^/x^ con- 
tribution in the B-term. Thus, they become relatively 
important at /i ~ 1. Another point is that as these terms 
arise as a non-linear coupling between the density and 
the velocity fields, they are of the order C)({i^']"(fc)}2). 
We thus may able to omit them on linear scales. 



B. Decomposition strategy 

We now describe how we decompose the power spec- 
trum in redshift space into spectra of density and velocity. 
Before that, let us introduce two useful quantities con- 
trolling the amplitude of the linear power spectra. We 
define 



P^^(fc,z) = 



--gUk,z)Pi^{k,zi..), 
gl{k,z)P^l{k,ziss), 



(8) 
(9) 



where ziss stands for the redshift at the time of the 
last scattering. The parameters, gs and gs, describe 
the growth rate of the density and the velocity fields 
from that epoch. Since the spectrum Pg™{k,z\ss) is well 
constrained by observations of the CMB temperature 
anisotropy, the parameters solely capture the properties 
of the gravitational law, and are expected to be free from 
the assumptions in the initial condition. 

In reconstructing the spectra, we first bin the power 
spectrum measured from simulations into bins of k and fi. 
Then, for the i-th bin of wavcnumber fc, which we denote 
ki, we estimate Ps^s^ih), PgQ{ki) and PQQ{ki) based on 
the fj, dependence. We assume that linear 6 and linear Q 



PmPbb^ and treat two 



are perfectly correlated, Pjg ~ \/ 55 ee ' 
parameters, gsiki) and gQ{ki), as free parameters for the 
fci-bin. In fitting P{ki,^), we scale the spectra according 
to the parameters gs and g^: 
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fid.lii 



(fc). 
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where quantities with subscript "fid" are computed for 
the fiducial cosmological model. 

When we restrict the analysis to linear regime, the 
above procedure is expected to work properly. We further 
elaborate the procedure to correctly handle the nonlin- 
earity. We adopt the closure approximation (Taruya & 
Hiramatsu (3l| ) to predict the nonlinear growth in the 



three spectra, Pijik). In doing so, we simply assume the 
fiducial GR cosmology used in running the simulations. 
We pre-compute the nonlinear corrections to the three 
spectra up to the second-order in the Born approxima- 
tion: 



5Pt^{k)^Pt^{k)-P^^''\k), (11) 

while we allow to vary the linear part according to 
Eq. (jlOp . By adding up linear and nonlinear parts, we 
have 

P.i^) = ll) |^|^''"(^) + (12) 

Strictly speaking, the second term depends on the cos- 
mological model as well as the gravitational law. As a 
first trial, however, we simply let this term unchanged 
from its fiducial value. In Sec. IIV[ we will relax the as- 
sumption to execute a more general analysis: we adopt a 
wrong cosmological model as the fiducial model, and see 
how well we can recover the true spectra. 

We also scale the correction terms, A{k, fi) and B{k, ji) 
as follows. In every step of the fitting, given set of gs{k) 
and gQ{k), we compute their simple arithmetic means: 

95^ i^^9s{ki), 9e = i^^9e{ki), (13) 

■i I 

where the subscript i runs over fc-bins up to a maximum 
wavenumber, fcmax, and we denote the number of the bins 
by A^bin- Again, we pre-compute the correction terms for 
the fiducial cosmology using the standard perturbation 
theory, which we denote A^'^ and B^'^, and scale them 
according to the average values of gs and gg : 



A{k,fi) = A 



fid 



B{k,fi) = B 



fid 



, 9s 99 

;;fid' ;;fid 

9s 9e 
, 9s 96 



(14) 



In the above, the terms originated from 6 (8) are multi- 
pUed by gs/gl'' {ge/gl")- 

We finally explain our strategy for the FoG factor. We 
try both Gaussian and Lorentzian functions, and we let 
ay as a free parameter. This parameter is determined by 
fitting globally the broadband shape of P{k,fi). 

In summary, our reconstruction strategy is as follows. 
We model the redshift-space power spectrum, P{k, fi), as 



P{k,fi) = 




9s{k) 
9lHk) 



gs{k) gejk) fj^.n, 
9r{k)9%Hk) 
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FIG. 1: The observed spectra of P{k,fi) are presented at scales 
0.03/iMpc~^ < fc < 0.11 hMpc^^. Blue triangles represent av- 
eraged P(fe, fi) at given k-fi bins from 60 realizations. The error 
bars are dispersions of measured values of 60 realizations divided 
by ^/60. 



fast Fourier transformation (FFT) using the cloud-in- 
cell (CIC) method. We use bins in k and ^ for the 
fohowing analysis, k is divided in Afc = 0.01/iMpc~^ 
hnearly equally spaced bins from k = 0.01 /iMpc"^ to 
0.2/iMpc~^ and fj. is in 20 linear-bins from to 1 with 
equal spacing. The averages of measured 2D power spec- 
tra in (/c./i) coordinate are shown in Fig.[TJ The Gaussian 
variance is used to derive errors for each bin shown as er- 
ror bars in Fig. [D a[Pob{k,fi)] = Pob{k, ii)y^2/N{k, fi) 
where Nik^n) is number of modes in 2048^( /i^^ Mpc)^ 
in Fourier space. 

The overall amplitude of P{ki,^j) at /i — > is solely 
determined by Pgg^{ki). And the running of P{ki,Hj) 
along ^ direction is determined by PQ^{ki) about the 
pivot point of P{ki,^j ~ 0). Those distinct contribu- 
tions of Pg™ and Pqq to the observed spectra lead us to 
simultaneously decompose both through data fitting to 
the observed P{ki,^j) in k and ^ dimension. Addition- 
ally, non-perturbative effect is externally parameterized 
using (T^ appearing in , while the higher order loop 
corrections are expressed using the given linear spectra 
parameters of P]™ and Pqq- Here parameterised one- 
dimensional velocity dispersion is set to be a scale- 
independent free parameter. 

We find best-fit parameter space of Pgg{ki), PQQ{ki) 
and ay by minimizing p9| . 



where factors with the subscript "fid" are pre-computed 
for the fiducial cosmology, and gs{ki) and ge{ki) as well 
as (T„ are free parameters. We show results for various 
cases in what follows. We include/exclude A{k,^) and 
B(k,fi), we adopt Gaussian and Lorentzian for the FoG 
damping factor, and we vary the maximum wavenumber 
included in the analysis, fcmax- 



III. MEASUREMENTS OF LINEAR SPECTRA 

A. The observed spectra P{k, from mock 
catalogues 



i—i^^i^ p—1 q—1 

X CoVp^^(fc,)[Pob(A:i,/ig) - P&tih.lJ-q)] , (16) 

where fcmin is fixed to be fcmin — 0.01 /iMpc^^, and best 
fcmax is determined in the following subsections. Off di- 
agonal elements of the covariance matrix are nearly neg- 
ligible and those diagonal elements are written as 



Gov p (fc,) 



1 



(17) 



We repeat this procedure for each 60 realization, and 
report the averages of best-fit values. 



For subsequent analysis, we use the dark matter distri- 
butions created by the simulations in Ref. [32|. The vol- 
ume size of the A''-body simulations is (2048ft.~^Mpc)^, 
and we have 60 independent snapshots at each of the 
four output redshifts, z = 0.35, 1, 2 and 3. The fiducial 
cosmological parameters of the simulation are given by 
(fl„ = 0.279, flf, = 0.165 17™, r^fc = 0,h = 0.701,0-8 = 
0.816, TT-s = 0.96). The distribution of dark matter par- 
ticles is modified according to their peculiar velocity to 
incorporate the redshift distortion effect. We adopt the 
distant-observer approximation and measure the power 
spectrum in {k±,k\\) space, where subscripts '_L' and 
'II' denote perpendicular and parallel components to the 
line-of-sight. 

The density fluctuation field is constructed by assign- 
ing the dark matter particles to 1024^ grids for the 



B. Decomposition of linear density— density spectra 

We present the result of decomposition of Pg^g^- Lin- 
ear spectra of Pjjj" are determined at a few first bins 
about /i = in which the orientation of correlated two 
point pairs is transverse. The observed Pob at those fi 
bins are nearly equivalent to density-density spectra it- 
self. The procedure to decompose P]™ is immune from 
all line-of-sight contamination described in SecHIl There- 
fore density-density spectra are measured in high preci- 
sion at arbitrary scale of k. But what we are interested 
in is measuring linear spectra of P]'/. Unless the non- 
linear contribution in the decomposed Pgg are separated, 
cross-correlation between density and velocity fields is 
not guaranteed to be perfect. 
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FIG. 2: Decomposed are presented with fcmax = 0.2hMpc ^. 
(upper panel) Dashed curve represents fiducial Pj'^"- Blue triangles 
represent fitting results using theoretical formulation in Eq. 1151 and 
black circles represent fitting results without non-linear correction 
terms in Eq. 1151 [bottom panel) The fractional errors are presented 
for measured Pg^^ in the upper panel. 



FIG. 3: Decomposed Pqq are presented with varying fcmax- (upper 
panel) Solid curve represents fiducial P@@. Blue triangles represent 
fitting results using fcmax = 0.11 hMpc~^, and black circles repre- 
sent fitting results using fcmax = 0.15/iMpc~^. (bottom panel) 
The fractional errors are presented for measured Pqq in the upper 
panel. Thin blue curves represent the estimated errors using the 
Fisher matrix analysis. 



The closure approximation in Eq. is applied for 
extracting linear information out of measured Pss- This 
approximation breaks down at specific scale of k. The 
upper bound of k is investigated to uncover the limit of 
the closure approximation for density fields. 

We use measured Poh{k,fJ,) up to k = 0.2/iMpc~^ at 
z = 0.35. Then there are 400 measured Poh{k,ij) at each 
realization of simulated maps. Those measured data are 
fitted by the parameter space of {2QPf^, 20P^5^, a^). 
In order to handle a large parameter space, we adapt 
Markov chain Monte Carlo (MCMC) methods which are 
a class of algorithms for sampling from probability dis- 
tributions around equilibrium points. 

Dashed curve in Fig. [2] represents the fiducial P]'^"(/c) 
in logarithmic scale. Black circle points in the top panel 
represent measured Pss{ki) when non-hnear correction 
terms of 5Pij in Eq. (jlip are nullified. Fractional errors 
against linear spectra of density fields are shown in the 
bottom panel of Fig. [2] Non-linearity of density fields 
in measured Pss{ki) is observed from k = 0.1/iMpc~^. 
Measured P]^"(fc.O deviates from Pffik) by 10% at k = 
0.2ft,Mpc~^. Blue triangle points represent measured 
Pss{ki) including 5Pij. Seen at fractional error bars 
in the bottom panel, linear spectra of are well 

reproduced within a couple of percentage uncertainties 
through k ~ 0.2ft.Mpc~^. The closure approximation in 
Eq. pT|) for density fields is proved to be trustable at 
least by this limit of k. 

Cross-correlation coefficient between density and ve- 



locity fields is not parameterized in this paper. Instead, 
we claim to probe linear spectra of density and veloc- 
ity fields using this closure approximation. Our test for 
linearity of measured density spectra in this subsection 
is important to proceed next step of probing coherent 
motion spectra. 

Additionally, measuring Pgg^{k) is also precious infor- 
mation to probe distance measures. The important infor- 
mation about the evolution of the universe is imprinted 
on the large-scale structure of the universe. The broad- 
band shape of the power spectrum provides the informa- 
tion about the horizon scale at the epoch of the matter- 
radiation equality [ssl - fSTj . Although much of this infor- 
mation has been faded away due to the nonlinear process, 
the possible extension of distance measure detectability 
at bigger k is worth being investigated in the following 
works. 



C. Decomposition of coherent motion spectra 

Results of Pqq decomposition are presented in this 
subsection. The running of Pob alon_g the fi direction is 
caused by Pee' ^^'^ pivoted from Pob at /i = 0. The 
observed Pob are maximally affected by peculiar velocity 
when the orientation of correlation is radial, while mini- 
mally affected at transverse orientation of correlation. In 
the limit of ^ — 1, Pob is significantly contaminated by 
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FIG. 4: Correlation between Pqq and av is presented at the bin 
of A; = 0.10/iMpc~^, using one realization of simulated catalogue 
out of 60. The upper bound of k = O.ll/iMpc-i. It 

corresponds to the last bin of blue triangle points in Fig. [S] 



all non-linear smearing effects described in SecHH Hence, 
we develop an appropriate statistical tool to bridge the 
improved theoretical model to real datascts. 

The statistical methodology to treat non-linear cor- 
rection terms is described in Sec. IIIBI in detail. In this 
method, the higher order loop correction terms are given 
theoretically, but the uncertainty due to FoG effect is 
phenomenologically parameterized. At scales in which 
the first order term of FoG dominates, the assumption of 
FoG effect is valid through fitting with datasets. But 
beyond this quasi-linear cut-off, our statistical model is 
broken down. The decomposition of Pqq in our analy- 
sis is limited by the uncertainty of FoG effect at higher 
orders. 

Using the simulated maps at z = 0.35, we find this 
upper bound of k scales in which our assumption is 
valid. We test two different upper bounds of fc,„ax — 
0.11 /iMpc~^ and 0.15 /iMpc"\ The blue triangle points 
in Fig. [3] represent the decomposed Pqq using /cmax = 
0.11 ft,Mpc^^. The fiducial spectra are well decomposed. 
The best fit is 3.7 h^^ Mpc with fc^ax = 0.11 /i Mpc"^ 
which corresponds to about 10% non-perturbative con- 
tribution to Poh at /i = 1. The detailed functional form 
of FoG is not crucial in our statistical method, because 
most FoG functions agree at the first order approxima- 
tion. The decomposed P^Qih) is not much dependent on 
types of FoG function, such as Gaussian or Lorentzian, 
in this test. Hereafter, Gaussian function is chosen for 
describing FoG effect. 



FIG. 5: We test the contribution of higher order terms of A(k,fi) 
and B{k, fj.) in Eq. [15] The upper bounds of k is fixed at /cmax — 
O.llhMpc"^. Dashed curves represent the fiducial Pqq, and blue 
triangles are fitting results using theoretical formulation in Eq. 1151 
Black circles in the upper panel represents fitting results without 
bi-spectral higher order terms of A{k, fi) in Eq. 1151 black circles 
in the bottom panel represents fitting results without quadratic 
higher order terms of B(k,ii) in Eg. 1151 



the measured PQQ{k) becomes underestimated [sO]- The 
black circle points in Fig. [3] represents decomposed 
Peeik) with fc,nax = 0.15/iMpc^^ This test indicates 
that the unknown higher order terms of FoG effect be- 
come dominating above k = 0.11/iMpc~^. Beyond 
^max = 0.11/iMpc~^, the first order approximation of 
FoG effect is not valid, and the decomposed Pqq start to 
be biased. 

Wc compare fractional errors of decomposed PQQ(fc) 
with the theoretical estimation using Fisher matrix anal- 
ysis [l3l |. Wc do not marginalize the Fisher matrix 
with FoG effect. The estimated errors using fc^ax = 
0.11 /iMpc~^ are presented as thin blue curves in Fig. [S] 
When P0g|(fc) is not correlated much with FoG effect at 
smaller k ol k < 0.05/iMpc~^ , the observed and esti- 
mated errors agree to each other. But when PQQ(fc) is 
affected much by FoG effect aX k ^ 0.1 ZiMpc"^, the ob- 
served fractional errors increase. The correlation between 
Pee(fc = 0.10/iMpc-i) and ijy is presented in Fig 2] us- 
ing one realization of simulated maps out of 60. This high 
correlation causes the increasing fractional errors about 
a factor of 3. 

Results in Fig. [S] present the contribution of higher- 



When we increase fcmax beyond k^ 



0.11 /i Mpc" 



order polynomials in Eq. (|15p . The decomposed Pqq 
will be overestimated or underestimated without A or 
B. The fractional errors in the bottom panel of Fig. [5] 
show the results, when higher-order polynomials arc not 
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FIG. 6: Measured spectra of PQQ{ki) are presented at higher red- 
shift at 2 = 1, z = 2 and z = 3 in the top, middle and bottom 
panels respectively. 



considered at all. It is the case of common practice using 
the simple combination of Kaiser's model and FoG effect. 
We show that it is not sufficient for the future precision 
wide-deep surveys. 

Additionally, we test our method against higher red- 
shift data at z = 1 , 2 and 3 in Fig. [HI The upper bounds 
of fcniax for these redshifts are larger than the low red- 
shift of z = 0.35 as expected, since the damping scale of 
FoG effect becomes rather milder at higher redshift. We 
found that the linear velocity spectra Pqq are precisely 



measured up to /cmax — 0.12, 0.14 and 0.18 /iMpc 
z = 1, 2 and 3, respectively. 



IV. DISCUSSION 



at 



fully specify the cosmological model. This is indeed one 
of the main reasons why the precision measurements of 
the redshift-space distortions and/or baryon acoustic os- 
cillations are highly desired to unlock the nature of late- 
time cosmic acceleration. Hence, a large statistical error 
of the cosmological parameters would cause an erroneous 
estimation of the non-linear corrections that potentially 
leads to a biased decomposition of the coherent motion 
spectra. To quantify the size of this, we consider wrong 
cosmological models, in which the density parameter of 
the dark energy Qy differs from the fiducial value by 3, 
5 and 10%, keeping the equation-of-state parameter of 
dark energy fixed. To be precise, we adopt slightly larger 
values of ri„, while we fix the spatial curvature fife, the 
spectral index Us, the normalization of fluctuation am- 
plitude at CMB scales As , and some combinations of the 
parameters, fly^h^ and ^hh^-, since these are expected to 
be tightly constrained by the CMB observations. Then, 
we compute the non-linear corrections to the redshift- 
space power spectrum, and repeat the same analysis as 
examined in previous section. 

Fig. [7] shows the decomposition results of the coherent 
motion spectra at z = 0.35. Here, we plot the fractional 
errors of the resultant spectra in the cases adopting larger 
values of fly. As anticipated, there appears a clear sys- 
tematic trend that as increasing fly , the coherent motion 
spectra tend to deviate from the fiducial spectrum. Nev- 
ertheless, apart from k ~ 0.1 ft,Mpc^^, the resultant size 
of the bias is rather small. This is because we are basi- 
cally looking at the scales where the contribution of the 
corrections terms is small, and even the 10% change of 
the cosmological parameter gives a little effect on the 
non-linear corrections. Also, a slight mismatch of the 
non-linear corrections can be partly absorbed into the 
FoG damping factor, which further reduces the impact 
of wrong prior assumptions. Since the uncertainty of the 
future constraint on fii, will not be as large as 10%, the 
results shown in Fig. [7] may be regarded as a very good 
news, suggesting that the coherent motion spectra would 
be decomposed successfully in a less biased manner. 



The decomposition results shown in previous section 
rely on an idealistic assumption that the nonlinear cor- 
rections to the redshift-space power spectra, 6Pij, A and 
B, are known a priori for the fiducial cosmology. Since 
these corrections must be computed for a given set of cos- 
mological parameters, they might be sensitive to the un- 
derlying cosmological model and the decomposed power 
spectra could be biased if we adopt the wrong cosmologi- 
cal priors. In this section, we discuss the impact of wrong 
cosmological assumption on the decomposition of power 
spectra. 

In principle, the on-going and upcoming CMB exper- 
iments will provide a way to precisely determine a set 
of cosmological parameters, from which we can compute 
the non-linear corrections. Because of the parameter de- 
generacies, however, the CMB observation alone cannot 



V. CONCLUSION 

We have presented an improved prescription to recon- 
struct coherent motion spectra from the matter power 
spectrum in redshift space, properly taking account of 
the non-linear effects of both the structure growth and 
redshift distortions. Statistical analysis is presented to 
bridge theoretical models of redshift distortion to real 
datasct. Based on the perturbation theory treatment, 
non-linear correction terms of higher-order polynomials 
and non-linear growth functions are appropriately in- 
cluded in the reconstruction analysis. Those contribu- 
tions are proved to be infiuential even at linear scales, 
and measurements of coherent motion spectra are misled 
without it. On the other hand, non-perturbative cor- 
rection term such as FoG effect is parameterized. Those 
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FIG. 7: Fractional errors of measured spectra of PQQ(ki) are pre- 
sented at z = 0.35 using diverse templates at different cosmological 
models. Results with templates of Qv larger by 3%, 5% and 10% 
than true f!„ are shown at top, middle and bottom panels respec- 
tively. 



all knowns and unknowns in our analysis successfully re- 
produce fiducial spectra of coherent motion up to some 
limited scales. Although it is still steps away to achieve 
complete and practical observational tools measuring co- 
herent motion spectra, results show that our analysis 
method righteous path to be developed to meet real 
world. 
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